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Abstract. The analysis of the recent WMAP source catalog shows that the vast majority of bright foreground 
extragalactic sources detected in CMB maps are Blazars. In this paper we calculate the contamination of CMB 
anisotropy maps by this type of flat-spectrum, strongly variable and polarized extragalactic radio sources using up- 
to-date results from recent deep multi-frequency surveys. From a careful archive search and from multi-frequency 
catalog cross-correlations we found that more than 50 known Blazars (or Blazar candidates) expected to be above 
the sensitivity limit of the BOOMERANG experiment are included in the 90/150 GHz BOOMERANG anisotropy 
maps, a factor ^ 15 larger than previously reported. Using a recent derivation of the Blazar radio LogN-LogS 
we show that Blazars, whose counts continue to grow steeply down to faint fluxes, can sensitively affect CMB 
fluctuation maps with even moderate resolution and sensitivity. We calculate specifically that these sources induce 
an average sky brightness of 0.2Jy/deg 2 , corresponding to an average temperature of ~ 3 — 5 fiK. Moreover, we 
find that the associated level of fluctuations is of the order of Ci^Biazar = 1.3 x 10~ 2 fiK 2 sr at 41 GHz and 
Ce^Biazar = 6.5 x 10~ 4 fiK 2 sr at 94 GHz. Taking into account both Blazar variability, causing the detection of a 
number of weak sources that rise above the detector sensitivity during flares in long-exposure satellite experiments, 
and the many steep-spectrum radio sources that flatten at high frequencies, as well as the contribution of radio- 
galaxies, we find that the level of residual fluctuation due to discrete extragalactic foreground sources could be 
factor of ~ 2 — 3 higher than the above estimate. We show that the Blazar induced fluctuations contaminate the 
CMB spectrum at the level of » 20 - 50% at £ = 500 and 50 - 100% at £ = 800, and thus they cannot be neglected 
in the derivation of the primordial CMB fluctuation spectrum. Careful cleaning for Blazar contamination of high 
sensitivity /high resolution CMB maps is therefore necessary before firm conclusions about weak features, like 
secondary high-i! peaks of the CMB power spectrum or very weak signals like CMB polarization measurements, 
can be achieved. 
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1. Introduction 

The Cosmic Microwave Background (hereafter CMB) contains - in nuce - a large amount of cosmological information. 
The paradigm follows from the fact that the CMB power spectrum is sensitive to the cosmological parameters which 
describe the fundamental properties of the universe (see, e.g., Jungman et al. 1996, Bond et al. 1997). The ability 
of the CMB data to constrain the cosmological parameters hinges upon the existence of acoustic oscillations in the 
primordial plasma at the last scattering epoch and on the ability of CMB experiments to disentangle such primordial 
fluctuations from those induced by cosmic structures at lower redshift. 

CMB temperature fluctuations have been measured and studied in detail for the past decade starting with the 
large angular scale (0^9 deg corresponding to a scale I 20) COBE detection (Smooth et al. 1992) up to the more 
recent and small-scale detection by BOOMERANG (deBernardis et al. 2000), MAXIMA (Hanany et al. 2000), DASI 
(Halverson et al. 2002), VSA (Grainge et al. 2002), CBI (Pearson et al. 2002) and, lastly, by WMAP (Bennet et al. 2003) 
which probes the CMB angular spectrum till t ~ 800 over the whole sky with unprecedented precision. Extracting 
constraints on the cosmological parameters from the measurements of the CMB power spectrum is conceptually simple, 
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but much more delicate in practice. In fact, the ability of the present and upcoming CMB experiments to determine the 
cosmological parameters requires a careful cleaning of the CMB maps from the galactic and extra-galactic foregrounds 
(see, e.g., Refregier & Spergel 2000). 

At high Galactic latitudes point-like sources are one of the main contaminants at small angular scales where, 
even when the brightest sources are identified and then subtracted, their residual contribution may dominate over 
the cosmic variance uncertainty. It has been shown, in fact, that point-like sources are one of the main foreground 
issue for the CBI experiment which is sensitive up to £ ~ 3500 (Myers et al. 2002). However, such an issue also 
applies to other CMB experiments like VSA, DASI, or WMAP which probe the CMB angular spectrum at lower 
£ i5 1000 and it is nonetheless crucial for the up-coming space missions like PLANCK. Along this line of reasoning, 
Pierpaoli (2003) has recently shown that point-like sources induce a contamination of ~ 20% level in the CMB power 
spectrum reconstruction even at £ as 200 in the simulated WMAP map at 90 GHz. This level of contamination grows 
for increasing values of I (i.e., for smaller angular scales) reaching a level of ~ 80% at £ as 300. This implies that 
point-like sources must be successfully identified and efficiently subtracted off the CMB maps before recovering the 
true primordial fluctuation spectrum. Pierpaoli (2003) also showed that bright sources have the greatest impact on 
the power spectrum at moderate £ values {£ ~ 400) while the fainter sources have a greater impact at high-^ values. 

The recent WMAP experiment (Bennet et al. 2003) resolved 208 bright pointlike sources, the vast majority of which 
arc Blazars. We show in this paper that recent Blazar surveys, combined with detailed multi-frequency Spectral Energy 
Distribution analysis, provide evidence that a much larger number of Blazars and active galaxies must be present as 
positive temperature fluctuations already visible in the available CMB maps with even moderate spectral and spatial 
sensitivity. To this end, we discuss in Sect. 2 the astrophysical aspect of Blazars at radio and sub- mm frequencies 
and we discuss in Sect. 3 their contribution to the BOOMERANG 90/150 GHz anisotropy map. We also compare our 
results with the catalog of bright sources derived from the WMAP team. Based on an updated derivation of the Blazar 
counts, we derive in Sect. 4 the average sky brightness of the mm sky due to the Blazar population and we compute 
in Sect. 5 their angular power spectrum compared to the CMB power spectrum. We summarize our conclusions and 
we discuss the relevance of our results for the present and for the next coming CMB precision experiment in the final 
Sect. 6. 



2. Blazars at cm and mm wavelengths 

Blazars are a type of radio loud Active Galactic Nuclei (AGN) distinguished by their extreme, in some cases unique, 
properties. These include large amplitude rapid variability, apparent superluminal motion, high level of polarization 
and emission of radiation across the entire electromagnetic (e.m.) spectrum, from low frequency radio waves to gamma 
rays, up to TeV energies in some cases. This peculiar behaviour is thought to be due to physical phenomena observed 
under very special circumstances such as the emission of radiation from a relativistic jet of material ejected by a radio 
galaxy in the direction of the observer and seen at small angle with respect to the line of sight jUrry fc Padovani 19951 
IBlandford fc Rees 19780 . 

The broad band e.m. spectrum of Blazars [represented in a vf{v) vs. v plane and usually referred to as the Spectral 
Energy Distribution or SED] is marked by a flat radio spectrum, widely believed to be due to synchrotron radiation, 
whose luminosity increases with frequency until it reaches a maximum that is normally located at 10 12 — 10 Hz 
but can also reach ~ 10 17 Hz in some extreme cases. Above the synchrotron peak the spectrum falls sharply until 
Inverse Compton (IC) emission emerges and rises again to form a second luminosity peak which is typically located 
at frequencies v as 10 17 — 10 19 Hz [see, e.g., FigQlfbr two examples of Blazar SEDs; many other examples of Blazar 
SEDs can be found in the catalog of Giommi et al. (2002a)| . 



Because the special conditions that make a radio galaxy appear like a Blazar are rather improbable, Blazars are 
much rarer than other AGN. However, since Blazars emit over the entire e.m. spectrum, while other types of AGN 
emit most of their radiation in the radio (radio loud only), infra-red, optical, and X-ray band, Blazars remain the 
extragalactic sources more frequently found in the remaining parts of the e.m. spectrum, that is the Gamma-Ray, the 
TeV and the microwave bands. 

There are two types of Blazars: those showing a featureless spectrum, called BL Lacs, which are rare and do 
not show a strong cosmological evolution, and the flat radio spectrum QSOs (FSRQ), which have broad lines in 
their optical spectrum, are much more abundant and evolve strongly on cosmological times | |Urry fc Padovani 19951 
IPadovani 20020 . Intensity and spectral variations are conspicuous at optical frequencies and become even larger at 
higher frequencies. Variability in the radio part of the e.m. spectrum, as demonstrated by a monitoring program lasting 
more than 25 years, can reach factors of ~ 2 to 10 on typical time-scales of weeks to years with a clear tendency to 
become more pronounced and more frequent at higher frequencies; polarization at radio frequencies is often present at 
a level of a few percent and in some case up to 10% (see, e.g., [Hughes et al. 1992| and the on-line UMRAO database 
http: / /www. astro. lsa.umich.edu/obs / radiotel/umrao.html) . 
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Fig. 1. Left We show the SED of the FSRQ source, PKS 0438-436, one of the three Blazars reported by the original 



authors to have been detected in BOOMERANG maps (deBernardis et al. 2000, Coble et al. (2003) I. The three open 
diamonds at 90, 150 and 240 GHz arc, in fact, BOOMERANG measurements taken from |Coble et al. (2003)| The 



dotted and dashed-dotted curves below the data points, represent the frequency spectrum of CMB fluctuations of 
the order of 300 and 50 [iK , respectively. The solid curve shows the SED in the v — /„ plane as expected from the 
Synchrotron Self Compton (SSC) model. Right We show here the SED of the FSRQ source PKS 0454-463, whose 
expected flux at mm frequencies, estimated from the SSC model and shown as a dashed curve, is well above the 
BOOMERANG detection limit. This extrapolation is confirmed by the WMAP detection of this source at 23, 33, 41, 
61 and 94 GHz ^Bennett et al. 2003j) which are shown as open points. The dotted and dashed-dotted curves below the 
data points, represent again the frequency spectrum of CMB fluctuations of the order of 300 and 50 [iK, respectively. 



Radio sources of this type display a spectral shape which is not too different from that of the CMB spectrum before 
and around the peak of its power in the frequency range ~ 50 — 200 GHz (see, e.g., Fig.l). For this reason Blazars 
can easily contaminate CMB maps (e.g., they can be confused with CMB fluctuations) if they are bright enough to 
be detected by CMB experiments. The recent results of the WMAP satellite experiment (jBennett et al. 2003|l clearly 
demonstrate that Blazars are present in large numbers in CMB fluctuation maps. In fact, we found, from a cross- 
correlation of the 208 sources of the WMAP catalog with catalogs of known sources and with the NED database, 
that ~ 85% of the WMAP sources are Blazar-like objects. Specifically, we found 139 FSRQ, 24 BL Lacs, and 15 radio 
galaxies. In addition to this evidence, we found that the three extragalactic sources detected well above the sensitivity 
threshold of the BOOMERANG maps IjDe Bernardis et al. 2000l|Coble et al. (2003)|) are also Blazars. 



To produce a more quantitative analysis of the contamination of CMB anisotropy maps by Blazars, we provide 
in the following, a list of the bright known Blazars which can be found in the 90/150 GHz BOOMERANG map. 
We choose to analyze the BOOMERANG maps because it covers a statistically representative large area of the 
sky ~ 1800 deg 2 ) and because they are available (even though in a low resolution format) on the web (see, e.g., 
http: / /cmb. phys.cwru.edu/boomerang/ ) . 

3. Blazars in the BOOMERANG CMB anisotropy maps 

We search here for Blazars that might have already been detected in the BOOMERANG maps but not yet recognized 
as such. To this aim, we have built the SED of all known Blazars for which that are enough published measurements 
(taken from NED and from various other surveys) to allow for a reasonably safe extrapolation of the flux intensity at 
the frequencies 90 and 150 GHz of the BOOMERANG maps. We started our analysis from the few point-like sources 
that have already been recognized in the BOOMERANG maps l|De Bernardis et al. 2000)) . Their follow-up (Coble et al. 
2003) has shown their galactic (6 out of 9 sources) or extra-galactic (3 out of 9 sources) nature. The three extragalactic 
foreground sources are namely a Blazar (PKS0537-441), a BL Lac (PKS0521-365) and a Flat Spectrum Radio QSO 
(PKS0438-43). For one of these, namely PKS0438-436, we report the SED in FigH where filled circles represent NED 
measurements and the three open diamonds at 90, 150 and 240 GHz represent the BOOMERANG measurements 



taken from Coble et al. (2003) To obtain a flux estimate around the peak of the CMB we fitted the multi-frequency 
points with a Synchrotron Self Compton (SSC) model ( |Perri et al. (2002)) see also Perri 2003, PhD Thesis, University 
of Rome, La Sapienza) which is a very good approximation of the SED of Blazars in the radio and mm band. The 
dotted and dashed-dotted curves in Fig represent the spectral distribution of typical CMB temperature fluctuations 
of the order of 300 and 50 [iK respectively, which can be clearly detectable in the BOOMERANG maps. Based on 
the SED shown in Fig^ the source PKS0438-436 is then associated with a fluctuation of order ^ 300 \xK, as already 
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recognized in the BOOMERANG maps. The SED of the other sources PKS 0521-365 and PKS 0537-441 are similar 
to the previous case but are not reported here for the sake of brevity. 

Following a similar procedure we found more than 50 Blazars, or Blazar candidates, selected through a multi- 
frequency technique like that used in the DXRBS IjLandt et al. 20011 IFaTfovani 20021 using the NVSS (radio), GSC2 
(optical) and the Rosat (X-ray) all sky survey, for which the extrapolation of the radio data is above the detection limit 
of the BOOMERANG experiment. Since the SED of Blazars allows for fluxes above the BOOMERANG sensitivity 
and below the Rosat all sky survey sensitivity limit, the number of objects found is only a lower limit. 

As an example of our findings, we show in the right panel of Fig^ the SED of the source PKS 0454-483 whose 
expected emission at mm frequencies is similar to that of the three Blazars already found in the BOOMERANG map. 
Whenever a Blazar was found to have a flux corresponding to a CMB fluctuation significantly above 50 fj,K, we have 
identified its position in the 150 GHz BOOMERANG map. 

We report in FigHthe location of the whole set of Blazars found in the 150 GHz BOOMERANG map. Their sky 
coordinates (RA and DEC), radio flux, optical V-band magnitude and redshift are also reported for completeness in 
TabQ] To summarize, we found 54 known Blazars with flux greater than ~ 200 mJy that are probably associated 
with fluctuations of amplitude ft 50/iif, in the 150 GHz BOOMERANG map. The number of these Blazar point-like 
sources is a factor > 15 larger than previously reported ( |Coble et al. (2003) I. A much larger number of fainter Blazars 
are evidently expected to be associated with lower amplitude ( ^ 50/zif) CMB fluctuations due to the fact that the 
Blazar counts steeply increase at flux fainter than 200 mJy, as will be discussed in the next Sect. 4. 

The analysis of the BOOMERANG map is just an example of what we could find in a specific area of the sky sur- 
veyed at mm wavelenghts. In particular, we estimate that more 200 Blazars should be found above the BOOMERANG 
sensitivity limit in 1000 deg 2 area. The WMAP anisotropy maps have basically the same sensitivity threshold and 
spatial resolution of the BOOMERANG maps and thus our results can be quite easily rescaled to the WMAP all-sky 
survey. Given the previous results, the Blazar populations are expected to produce a non-negligible contamination of 
the CMB maps as well as of the CMB power spectrum. We will address this point more specifically in the next Sects. 4 
and 5. 



4. The sky brightness due to Blazars at cm and mm wavelengths 

Precise estimates of the contribution of Blazars to the sky brightness at radio and mm frequencies have so far been 
hampered by the lack of sufficiently deep surveys. For a long time only very shallow (1 or 2 Jy) surveys have been 
available at cm frequencies IjStickel et al. 199"TllWall fc Peacock 1985|l . However, the recent multi-frequency approach 
to Blazar search is changing the picture (e.g. ILaurent^M uchleis en et al. 19981 IGiommi et al. 19991 [Padovani 20021 
IGiommi et al. 2002blGiommi et al. 2002cjl . making it possible to detect many fainter Blazars. 

The integrated radio flux due to a population of discrete sources in the flux range S ra j„ — S max can be calculated 
as follows: 

where the quantity dN/dS is usually referred to as the differential LogN-LogS distribution that, for a population of 
sources with constant density, is given by a power-law with slope equal to 2.5 which then flattens at faint fluxes of 
order of ^1 — 10 mJy, depending on the Blazar luminosity function. The LogN-LogS of BL Lacs is presently known 
down to ~ 3.5 mJy for the subsample of extreme High Energy Synchrotron peaked (HBL) sources IjGiommi et al. 19991 
IGiommi et al. 2002b) while the best estimate of the LogN-LogS of the full population of BL Lacs and FSRQs currently 
comes from the Deep X-Ray Blazar Survey (DXRBS. ILandt et al. 200l)l . This is the deepest complete and radio-flux 
limited Blazar survey presently available which extends to 50 mJy (Pad ovani et al. 2003(1 . At radio frequencies, FSRQs 
are more abundant than BL Lacs by about a factor of « 4 and show a strong cosmological evolution, while BL Lacs 
cither do not evolve or show negative cosmological evolution iPadovani 20021 IGiommi et al. 19991 IBade et al. 19981 
IRector et al. 2007)1 1 . The LogN-LogS of BL Lacs is approximately euclidean up to at least 50 mJy, while below 20 mJy 
it flattens substantially, at least in the subsample of extreme HBL BL Lacs (Gio mmi et al. 2002b|l . From the DXRBS 
survey we know that the surface density of sources brighter than 0.1 Jy at 5 GHz is 0.1 object per square degree for 
FSRQs and 0.025 objects per square degree for BL Lacs (Giomm i et al. 2002bMPadovani et al. 2003^1 . We can then 
assume a Blazar overall density of 0.125 deg" 2 with flux S > 0.1 Jy. Since FSRQs (the dominant Blazar population) 
evolve strongly, while BL Lacs show no evolution (or slightly positive) at high radio fluxes l|Padovani 2 002'). and 
negative evolution below 20 mJy IjOiommi et al. 19991 IGiommi et al. 2002bl) , we can conservatively assume an overall 
Blazar differential LogN-LogS with an "Euclidean" slope of 2.5 which extends up to a break flux of the order of 10 
mJy below which it flattens significantly to a slope of ~ 2.2. This assumption is fully supported by the LogN-logS 
(see FigOU that we have built using a sample of 135 Blazars (BL Lacs + FSRQs) from the NVSS-RASS Uy and the 
0.1 and 0.05 Jy points of the DXRBS l|Padovani et al. 2003|) . While this LogN-logS is consistent with the counts of 



P. Giommi & S. Colafrancesco: Blazars in CMB maps 



5 



blazar_boorn_l 5 Qfihz.fi if 



Fig. 2. The 150 GHz CMB fluctuation map (taken from http://cmb.phys.cwru.edu/boomerang/) obtained from 
BOOMERANG is here correlated with the positions of known and candidate high galactic latitude > 20) Blazars 
(and radio galaxies) whose expected flux at 150 GHz (estimated from lower frequency archive data measurements) 
is high enough to be above the BOOMERANG detection limit, and should appear as green/yellow or red areas (see 
Table 1. See also http://www.asdc.asi.it/boomerang/ for the full list of FSRQ SEDs). Some of these objects (12 out 
of the 54 Blazars we list in Tab.l) are also included in the WMAP source catalog. 



flat spectrum radio sources derived by Toffolatti et aL (1998) at w 1 Jy, it does not show any flattening until a rather 



faint radio flux is reached and not just below 1 Jy as in the predictions based on Toffolatti et al. (1998) 



Thus, given the previously discussed LogN-LogS, and setting consistently the lower integration limit to 0.1 mJy, 
Eq.JEJ yields a radio brightness of w 0.7 — 1.0 Jy/deg 2 at 5 GHz. In order to properly extrapolate this brightness to 
higher frequencies, we have to make use of the distribution of radio to mm spectral slopes for the Blazar population. 
To estimate realistic slope values we have built the Spectral Energy Distribution (SED) of all Blazars in the NVSS- 
RASS Uy survey IJCiommi et al. 2002c) using the many radio measurements available from the literature (mostly 
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Fig. 3. Left. The 5 GHz (integral) LogN-Logs of Blazars built with data from the NVSS-RASSS Uy survey and the 
DXRBS point at 0.1 and 0.05 Jy. The dashed line is a power-law with slope 1.6. Right. We show here the distribution 
of spectral slopes between cm and mm wavelenghts from the sample of 135 Blazars of the NVSS-RASS Uy survey 
with a < 0.5. 



taken from NED). Since all objects in the NVSS-RASS survey are brighter than 1 Jy, several radio data points are 
normally available for each Blazar over a wide frequency range. We have then fitted a a power law to the radio part 
of spectrum and measured the distribution of spectral slopes a [here f(v) oc v~ a \ between the cm (e.g., 5GHz) 
and mm region (e.g., « 41,90,94 and 150 GHz). We have derived the normalized distribution shown in Fig|3] using 
the subsample of 135 objects with a < 0.5 which is a requirement normally applied to build Blazar samples. Using 
such a spectral slope distribution, the expected 150 GHz sky brightness due to Blazars ranges from « 7 x 10 10 to 
« 1 x 10 11 Jy Hz deg~ 2 depending mildly on the assumptions for the LogN-LogS parameters. Since the CMB brightness 
at 150 GHz is « 2.1 10 16 Jy Hz deg~ 2 (e.g. |Hasinger 2000| ), the Blazars brightness is w (3 - 5) x 10~ 6 that of the 
CMB, which is only a factor of a few below the level of the presently measured CMB fluctuations and similar to the 
sensitivity planned for the upcoming PLANCK mission. At 41 and 94 GHz we obtain a sky brightness of 3-5 fj,K 
assuming that the slope of the LogN-LogS flattens from 2.5 to 2.1 at 10 mJy and integrating it down to 1 and 0.1 
mJy, respectively. 

This estimate is only a lower limit to the average sky brightness induced by Blazars since the DXRBS and NVSS- 
RASS Uy surveys only include sources with flat radio spectrum between 1.4 and 5 GHz, while it is well known that a 
number of radio source spectra that are steep at low frequencies progressively flatten at higher frequencies. Through 
direct 15 GHz source counts measurements (see Taylor et al. 2003) estimate that the contribution of steep-spectrum 
sources at 1.4-5 GHz increase the extrapolated source counts at high frequency by a factor of ~ 2. 

In addition, other sources that are not accounted for in Blazar surveys have been detected, as a minority population 
but still in fair numbers, in the WMAP experiment. These are, e.g., radio galaxies which have a often extended emission 
with steep radio spectrum, that outshines a much flatter spectrum nuclear component which may emerge and dominate 
the emission above a few tens of GHz. One example of such sources is Pictor A, a bright radio galaxy with extended 
steep-spectrum emission and a flat-spectrum Blazar like, nuclear component (Fiocchi & Grandi 2003, in preparation), 
which is included in the WMAP source catalog and is also present in the BOOMERANG maps (see Sect. 2). 



5. The angular power-spectrum contributed by Blazars 

Assuming that the Blazar are spatially distributed like a Poissonian random sample, then the angular power spectrum, 
Ce, Blazar, contributed by these sources is given by 

Ge, Blazar = J dS S 2 , (2) 

(Tegmark and Efstathiou 1996, Scott and White 1999). The quantity at right hand side in Eq.(|2J) is the usual Poisson 
shot-noise term (Peebles 1980, Tegmark and Efstathiou 1996) and a further term, lu(I) 2 , adds to it if the clustering of 
sources is not negligible (Scott and White 1999). This last clustering term can be even a factor ~ 10 larger than the 
Poissonian term at large and intermediate angular scales, £ < 1000 (Scott and White 1999). 

For the Blazar population described by the LogN-LogS given in Fig|2|we found Ci^Biazar ~ 31Jy 2 sr~ 1 at 41 GHz 
and Ce, Blazar ~ 30Jy 2 sr~ 1 at 94 GHz. Our results can be translated into temperature units, as is usual in CMB 
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Fig. 4. Left. We show the contribution of Blazars to the CMB fluctuation spectrum in the WMAP Q channel at 41 
GHz as evaluated from the LogN-LogS given in Sect. 4 (solid line). We also show the angular power spectrum for the 
Blazar population by adding the contribution of radio sources with steep-spectrum at low radio frequencies which 
flatten at higher frequencies (dashed line). The dotted line includes also the effect of spectral and flux variability (see 
text for details). A typical CMB power spectrum evaluated in a ACDM cosmology with fl m = 0.3, = 0.65, fib = 0.05 
which best fits the available data is shown for comparison. Right. We show here the ratio Ci^Biazar/Cg^cMB between 
the Blazar power spectra and the CMB power spectrum. 



anisotropy studies, using the conversion between the isotropic black-body (Planckian) brightness B$(v) and the CMB 
temperature Tq which writes as 



8Bq k ( kTo 



dT 2 V he 



2 r 9. n 2 / nA n t -1\ r ,2 



X 



sinh(x/2) 



f 24.8 Jy sr' 1 



x 



ih{x/2) 



(3) 



where x = hv/kT® — v/56.84GHz is the a-dimensional frequency given in terms of the Planck constant h, of the 
speed of light c and of the Boltzmann constant k. For WMAP, we found Cg^Biazar ~ 1-3 • 10~ 2 [iK 2 sr at 41 GHz and 
Ct^Biazar ~ 6.5 • 10~ 4 /iK 2 sr at 94 GHz. In Fig@]we plot the quantity Cg.Biazar^ + l)/(27r) and we compare it to 
the CMB fluctuation power spectrum which best fits the available data. The solid line in Fig0| shows the contribution 
of Blazars to the CMB fluctuation spectrum in the WMAP Q channel at 41 GHz as evaluated from the LogN-LogS 
given in Eq^ which yields Ct Blazar ~ 1-3 • 10 -2 /j,K 2 sr. We stress here that the previous value is a definite lower 
limit for Ci^Biazar since it neglects the contribution of steep-spectrum sources at low frequencies which flatten at these 
frequencies (41 and 94 GHz) and the effect of flux variability. For this reason we also show a more realistic case for 
the Blazar power spectrum by adding the contribution of radio sources with steep-spectrum at low radio frequencies 
which flatten at higher frequencies (the dashed line in FigQJ. These sources likely provide a value of Ci t Biazar which is 
higher by a factor ~ 2 (see also discussion in Taylor et al. 2003). The effect of spectral and flux variability also allows 
for an increase of the fluctuation level Ci^Biazar because many Blazars which are below the sensitivity threshold of 
CMB experiments can be detected during flares. The Blazar flux variability at millimeter wavelengths may be very 
substantial (higher than the factors 3-10 on time scales of weeks to years seen at cm wavelengths) and could definitely 
increase the level of contamination of CMB maps when these are built over long integration periods. The dotted line 
in Fig^ shows the Blazar power spectrum when also this last contribution is taken into account. 

To better quantify the level of contamination of the CMB power spectrum induced by Blazars we show in the 
right panel of Fig^Jthe ratio Ce_Biazar/C'g_cMB where Ci^cmb has been evaluated in a ACDM cosmology with O m = 
0.3, fl\ = 0.65, il — 0.05 which best fits the available data. It is evident from our results that, even in the case of the 
lower limit to Cg^Biazar (solid curve in Fig0J), the population of flat spectrum radio sources produces a power spectrum 
which amounts to ~ 20% of the CMB power spectrum at I i=a 500 and reaches ~ 50% of the CMB power spectrum 
at i pa 600. In this case the Blazar power spectrum equals the CMB power spectrum at I ^ 10 3 . In a more realistic 
case, the addition of steep-spectrum radio sources at low-i^ which flatten at higher z/s makes the contamination of 
the primordial CMB spectrum more significant. In fact, in this last case, the Blazar power spectrum (dashed curve in 
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FigEJ reaches a level ~ 40% at I w 500 and nearly reaches equality at £ ~ 600. Inclusion of variability effects increases 
even more the contamination (dotted curve in FigQJ reaching a level of ~ 60% at £ w 500 and equalizing the CMB 
power spectrum at £ ~ 550. Here we have not considered the statistical uncertainties of the CMB power spectrum data. 
These results show that Blazar-like (flat-spectrum) radio sources introduce a non-negligible contamination in CMB 
maps and in the CMB power spectrum reconstruction, even in the lower-limit case of considering just the flat-spectrum 
sources at high frequencies. Notice, that the Blazar contamination of CMB maps decreases with increasing frequency 
and it becomes < a few % at £ < 500 for v > 90 GHz. 

We also stress that the previous calculations are performed neglecting the clustering term cu(I) 2 and thus they 
must be considered again as a lower limit to the realistic angular power spectrum contributed by the Blazars. Such a 
clustering term could produce even a factor ~ 10 amplification of the power spectrum at any £ < 10 3 (see, e.g., Scott 
and White 1999). We will consider more specifically this additional contribution elsewhere (Giommi & Colafrancesco 
in preparation). 

It stems from the previous analysis, that the low-£ part of the reconstructed CMB power spectrum is heavily 
contaminated by Blazars even at £ ~ 400, where the second and third acoustic peaks are expected in the best-fit 
cosmological model (i.e., a scenario with J7 TO w 0.3,51a ~ 0.65,51b w 0.05). Thus, the non-negligible Blazar contam- 
ination of the CMB anisotropy spectrum at £ <; 400 has evidently an impact on the derivation of the cosmological 
information which is contained in the CMB anisotropy spectrum. The second acoustic peak at £ m 500 is considered 
as an incontrovertible evidence of inflationary sound waves (see, e.g., Hu & White 1997) and its amplitude is mainly 
sensitive to the values of Vt^h 2 and fl m h 2 . Increasing il^h 2 has the effect of decreasing the amplitude of the second 
peak relative to the first and third peaks. A contamination of ~ 20 — 50% at I w 500 would tend to decrease the power 
of the CMB spectrum at the position where the second peak is expected in a A CDM scenario. This effect, by itself, 
would bias the quantities flbh 2 and/or fl m h 2 toward high values. The third acoustic peak (and in principle all the 
acoustic peaks) is mainly depending on the amount of Dark Matter in the universe and its amplitude decreases for 
increasing values of il m h 2 . The high level of contamination (~ 50 — 100%) of the CMB power spectrum at £ w 800 as 
well as the intrinsic uncertainty of the data render the determination of the third acoustic peak problematic. The first 
acoustic peak at £ f» 300 is contaminated at a level 10% and hence the constraint on Slo set by CMB anisotropy 
data are little affected by the Blazar contamination, unless a substantial clustering of these sources can be considered 
as realistic. The present data do not provide strong insights on this last point and we will study it in more details in 
a forthcoming paper (Giommi & Colafrancesco, in preparation). 

Finally, we stress that the Blazar emission is significantly polarized at the level of ~ 2 — 10 % at radio and 
optical frequencies. As a consequence of the non-negligible Blazar contamination of the CMB maps, the primordial 
polarization patterns could be heavily contaminated by these polarized sources. Hence, a significant contamination 
of CMB fluctuation maps and the non-negligible fraction of residual unresolved sky brightness will complicate the 
detection of primordial polarization fluctuations, as these are expected to be present only at a very low level, i.e. at a 
level ~ 20 — 50 lower than the primary CMB fluctuations (see, e.g., Hu & White 1997). 

6. Conclusions 

In this paper we studied the contamination of CMB fluctuation maps induced by flat-spectrum extragalactic radio 
sources which shine in the frequency bands probed by BOOMERANG and WMAP. Blazars have in fact a typical 
spectrum which is similar to the CMB spectrum in the range 40 - 150 GHz. We predicted the level of brightness and 
of the induced CMB anisotropy using the most updated version of the Blazar LogN-LogS. 

In view of the results presented in the previous sections, the main conclusion of this paper can be summarized as 
follows: 

— Blazars are the major point-like, bright source contaminant of CMB maps. 

— We expect more than 100 Blazars to be present in the BOOMERANG 90/150 GHz maps above the sensitivity 
limit. We already identified more than 50 of these sources on the basis of a detailed analysis of known and candidate 
Blazars. We have also identified a large number of Blazars (178 out of a total of 208 sources) among the sources of 
the WMAP bright source catalog. 

— The importance of Blazars as contaminants of CMB fluctuation maps stems from the fact that their LogN-LogS 
is steeper than that of flat spectrum radio sources previously used to estimate the contamination of CMB maps 
by extragalactic sources. The observed Blazar LogN-LogS at 5 GHz combined with the measured spectral slope 
distribution implies a LogN-LogS at higher frequencies which is steeper than Euclidean at least down to 50 mJy 
at all frequencies relevant for CMB experiments. We expect, in fact, ~ 0.4 sources per deg 2 at F > 50 mJy, quite 
independently of the frequency. 

— We expect an induced fluctuation amplitude of the CMB map of order of C^Blazar ~ 31Jy 2 /sr which corresponds 
to a temperature fluctuation of ~ 0.013/iAf 2 sr at 41 GHz. The corresponding Poissonian spectrum of the Blazar 
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population equals the CMB power spectrum at I w 1000, and contaminates it at levels ~ 20% at i » 500 and 
- 50% at i » 800. 

— The previous estimates of Ct Blazar are a firm lower limit to the level of CMB fluctuation induced by the Blazar 
population. In fact, the value of the Poissonian term Cg t siazar is likely to be higher by a factor at least ^ 2 — 3 
considering both the Blazar variability and the contribution of steep-spectrum sources at low frequencies, radio- 
galaxies and possibly starburst galaxies at lower fluxes. We also stress that allowance for an appropriate clustering 
term might enhance substantially the angular power spectrum of the point-like sources and hence increase the level 
of contamination of the CMB power spectrum. 

— It is important to stress that flux variability at millimeter wavelengths may be very substantial (higher than 
the factors 3-10 on time scales of weeks to months seen at cm wavelengths) and definitely increases the level of 
contamination of CMB maps when these are built over long integration periods. Moreover, the Blazar variability 
increases with increasing frequency so that such a bias is more relevant for high-frequency CMB detections. 

— As a consequence of the non-negligible Blazar contamination of the CMB maps, the primordial CMB polarization 
seems to be greatly affected by the presence of polarized Blazars. Thus, the significant contamination of CMB 
fluctuation maps and the non-negligible fraction of residual unresolved sky brightness will render difficult to detect 
primordial polarization fluctuations. 

The many high-sensitivity and high-resolution CMB experiments planned for the near future will definitely face 
with the problems highlighted in this paper. In fact, we predict at least 1300 sources/sr at 90 GHz which produce a 
flux excess corresponding to > 10 \iK in a pixel of 15 arcmin size. One proven way of finding Blazars down to faint 
radio fluxes (~ mJy) is to use highly efficient multi-frequency selection techniques, similarly to the DXRBS survey 
which has an efficiency of ~ 90% l|Landt et al. 200Tjl . by cross correlating deep radio surveys (e.g., NVSS, FIRST etc.) 
with deep X-Ray data over the largest possible area of sky. However, to reach flux levels lower than ~ 10 mJy, this 
approach requires an X-ray sensitivity of ~ 10~ erg cm _2 s _1 . Such a sensitivity is reached in XMM-Newton and 
Chandra X-ray images only for a very limited area of sky. In conclusion, it is crucial to determine the level of Blazar 
contamination in precision CMB experiments in order to separate the contaminating foreground from the primordial 
CMB fluctuations and to provide a reliable determination of the cosmological parameters. 
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Table 1. The list of flat-spectrum radio sources (mostly Blazars) found in the BOOMERANG maps. For each object 
we list the source name (col.l), the coordinates Ra (col. 2) and Dec (col. 3), the radio flux at 5 GHz in units of mJy 
(col. 4), the optical V-band magnitude (col. 5), the redshift taken from NED (col. 6). Some notes on the listed sources 
are added in col. 7 of the table. 
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